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SUMMARY 


A  data  set  whkh  contains  nearly  500  arrays  representing  measurements  of  sky  and  terrain  radiance  values  is 
described  and  illustrated.  The  measurements  were  made  using  radiometer  systems  mounted  on  a  C-130  aircraft  during  a 
series  of  European  flints  associated  with  the  NATO  program  OPAQUE,  Fenn  (1978)  and  Johnson  et  ai  (1979). 

Radiance  measurements  throughout  the  4ir  field  surrounding  the  aircraft  were  made  in  four  spectral  bands  having 
mean  wavelengths  of  478,  557,  664  and  765nm.  Spectral  sets  of  measurements  were  made  at  each  of  several  altitudes 
between  the  surface  and  approximately  6  km  AGL. 

Evaluation  of  the  data  illustrates  good  orientation  and  radiometric  accuracies,  except  for  stray  light  influences  as  the 
radiometer  field  of  view  swept  near  and  through  the  solar  disc.  Corrective  procedures  are  outlined  and  their  application  is 
undergoing  further  refinement. 

These  data  are  appropriate  for  use  in  the  development  of  operationally  useful  predictive  models  as  illustrated  in 
Hering  (1981).  The  radiance  data  in  coryunction  with  their  companion  scattering  coefficient  data,  Johnson  and  Gordon 
(1980)  are  readily  applicable  to  the  determination  of  slant  path  contrast  transmittances,  atmospheric  optical  depths,  aerosol 
directional  scattering  characteristics,  flux  divergences  and  their  attendant  determinations  of  turbid  atmospheric  single 
scattering  albedos. 
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AIRBORNE  MEASUREMENTS  OF  EUROPEAN 
SKY  AND  TERRAIN  RADIANCES 


Richard  W.  Johnson 


1.  INTRODUCTION 

As  the  ever  broadening  world  of  electro-optical 
search  and  detection  continues  to  impose  its  technical 
influence  upon  our  most  imaginative  operational 
scenarios,  it  behooves  us  to  be  assured  of  our  abilities  to 
fully  understand,  and  hopefully  to  accurately  model  the 
characteriaation  of  these  operational  environments.  It  is 
in  support  of  this  general  context  that  the  Visibility 
Laboratory  in  cooperation  with,  and  under  the  sponsorship 
of  the  Air  Force  Geophysics  Laboratory  has  maintained  an 
extensive  program  of  airborne  optical  and  meteorological 
measurements.  Many  of  these  data,  all  of  which  were 
directed  toward  enhancing  our  understanding  the  funda¬ 
mental  elements  of  contrast  transmittance  through  the  tro¬ 
posphere,  have  been  presented  in  several  earlier  technical 
reports,  among  them  being  Dunlley,  etai  (1976),  John¬ 
son,  et  al.  (1979),  Johnson  and  Gordon  (1980)  and  John¬ 
son  (1981). 

Whereas  the  most  recent  of  the  reports  referenced 
above  have  concentrated  upon  the  documentation  and 
modelling  of  the  vertical  structure  of  the  atmospheric 
volume  scattering  coefficient,  it  is  important  to  note  that 
there  is  a  substantial  and  valuable  complementary  data  set 
associated  with  each  of  these  flight  episodes  that  has  not 
been  made  readily  available  to  the  user  community.  This 
complementary  data  set  includes  angularly  precise  spectral 
measurements  of  the  4ir  radiance  field  surrounding  the 
aircraft  at  several  different  altitudes  throughout  the  zero  to 
six  kilometer  altitude  regime.  The  possession  of  these 
radiance  data  is  what  enables  one  to  characterize  a  broad 
variety  of  the  environmental  conditions  extant  during  the 
flight  episode  and  thus  lead  to  the  development  of  opera¬ 
tionally  useful  predictive  models.  The  radiance  data  in 
conjunction  with  their  companion  scattering  coefficient 
data  are  readily  applicable  to  the  determination  of  slant 
path  contrast  transmittances,  atmospheric  optical  depths, 
aerosol  directional  scattering  characteristics,  flux  diver¬ 
gences  and  their  attendant  determinations  of  turbid  atmo¬ 
sphere  single  scattering  albedos. 

The  presentation  of  selected  sets  of  these  4ir  radi¬ 
ance  distributions,  and  a  discussion  of  several  pertinent 
characteristics  of  the  data  and  its  application  are  the  sub- 
jecu  of  this  report. 


2.  DATA  COLLECTION  SYSTEM 

The  airborne  data  collection  system  that  was  used  to 
cdlect  the  measurements  presmted  and  discussed  in  this 
report  has  been  described  previously.  Each  of  the  data 
reports  listed  in  Appendix  C  discusses  the  general  charac¬ 
teristics  of  the  specific  optical-meteorolo^cal  sub-systems 
appropriate  to  that  report's  data.  Similarly,  the  following 
paragraphs  will  briefly  review  the  pertinent  characteristics 
of  the  Automatic  27r  Scanner  Assemblies,  the  primary 
instruments  related  to  the  data  presented  in  this  report. 

An  annotated  illustration  of  the  total  airborne 
instrument  system  as  it  appeared  during  a  1977  deploy¬ 
ment  to  Lorient,  France  is  shown  in  Fig.  2-1.  Attention 
should  be  drawn  to  the  Upper  and  Lower  Hemisphere 
Scanner  locations  and  the  corresponding  locations  of  their 
companion  camera  systems.  A  more  detailed  view  of  an 
Automatic  Scanner  Assembly,  two  of  which  are  mounted 
on  the  aircraft,  is  shown  in  Fig.  2-2.  In  this  view  the 
instrument  has  been  removed  from  the  aircraft  and  is 
attached  to  its  ground  storage  rack. 

The  scanner  assembly  illustrated  in  Fig.  2-2  is  essen¬ 
tially  a  servo-following  telescope  assembly  associated  with 
both  an  optical  filter  selector  and  a  photomultiplier  detec¬ 
tor.  The  mechanical  orientations  associated  with  these 
ba»c  sub-assemblies  are  illustrated  in  Fig.  2-3  whidi  is  an 
annotated  outline  drawing  of  the  composite  assembly. 

The  azimuth  and  elevation  drives  in  the  2ir  Scarmer 
system  are  servo  follower  devices  whose  power  and  con¬ 
trol  circuits  reside  inside  the  cargo  compartment  of  the  air¬ 
craft  where  their  performance  can  be  monitored  by  the 
flight  crew.  A  typM  installation  of  the  control  console  is 
shown  in  Fig.  2-4,  and  described  briefly  in  the  following 
paragraphs  which  are  abstracted  in  part  from  Duntley, 
etai.  (1970 and  1972). 

The  airborne  scanner  control  console  consists  of  six 
individual  panda  mounted  in  a  standard  19  inch  tack 
(Fig.  2-4),  plus  a  remote  start/stop  function  tocated  at  a 
nearby  crew  station  (not  illustn^).  The  component 
panda  are:  (1)  scanner  control  panda,  two  each,  (2)  pro¬ 
grammer  control  pand,  (3)  progrunmer  pand,  (4)  gyro 
control  pand,  and  (5)  generd  purpose  power  disiribution 
pand. 
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n«.  2-4.  Automatic  2^  scanner  control  console. 


Scanner  Control  Paneis.  As  mentioned  earlier,  both 
the  scanner  azimuth  and  elevation  drive  sub-assemblies 
are  relatively  simple  servo  follower  systems.  Thus,  for 
each  scanner  assembly  as  shown  in  Fig.  2-2,  there  is  an 
associated  scanner  control  panel  which  contains  the  indivi¬ 
dual  demodulator,  servo  transmitter  and  drive  power  cir¬ 
cuitry  Spare  panels  were  always  on-board  for  in-flight 
replacement  in  the  event  of  electrical  malfunction. 

Programmer  Control  Panel.  In  order  to  direct  a  2tr 
scanner's  field  of  v  in  any  given  direction,  transmitter 
instructions  may  oe  generated  in  either  of  two  ways: 
(1)  manually,  via  graduated  hand  operated  dials  on  the 
programmer  control  panel  face,  or  (2)  automatically,  via 
the  solid  state  control  circuits  located  in  the  programmer 
panel.  Selector  switches  for  choosing  either  of  these  two 
control  modes  are  located  on  the  programmer  control 
panel.  Transfer  of  the  automatic  mode  start/stop  control 
to  the  remote  panel  at  the  crew  station  console  is  accom¬ 
plished  via  the  programmer  control  panel,  as  is  synchroni¬ 
zation  of  the  scanner  search  program  with  the  in-flight 
data  logger. 

Programmer  Panel,  The  actual  automatic  pattern 
that  the  scanner  field  of  view  traverses  is  establish^  by  a 
solid  state  logic  control  array  located  in  the  programmer 
panel.  In  the  early  development  stages  of  the  automatic 
2ir  scanner  assembly,  there  was  a  strong  demand  for  a 
flexible,  high  speed  pattern  control  assembly.  The  most 
appropriate  technique  available  at  the  time  involved  the 


dev^opment  of  an  in-house  design  utilizing  a  standaidized 
logic  card  and  rack  assembly  available  from  Computer 
Control  Co.,  Inc. 

The  original  pattern  called  for  a  full  2ir  steradian 
seard),  plus  reset,  in  32  seconds,  Duntley,  era/.  (1970). 
The  azimuth  drive  was  constant  at  one  revolution  per 
second.  The  devation  drive  followed  a  variable  rate  ramp 
functitm.  This  scan  rate  was  too  fast  for  many  applications 
and  consequmtly  a  timing  modification  was  inserted  to 
slow  the  entire  routine  down  to  one  fifth  its  original  speed. 
In  this  slower  mode  the  full  2ir  steradian  seardi  was 
accomplished  in  160  seconds,  with  the  azimuth  sweep  con¬ 
stant  at  5  seconds  per  revolution. 

With  the  advent  of  several  data  processing  sophisti¬ 
cations,  it  became  apparent  that  further  scan  pattern 
alterations  were  necessary.  Consequently,  in  early  1972  a 
nuyor  revision  to  programmer  panels  was  instituted.  The 
salient  features  of  this  modification  were  as  follows: 

A.  In  contrast  with  the  scheme  originally  established,  as 
in  Duntley,  etal.  (1970),  both  airborne  and 
ground-based  programmers  were  converted  to  gen¬ 
erate  the  same  seardi  pattern. 

B.  The  primary  azimuth  drives  were  slowed  to  10 
seconds  per  revolution.  On  the  airborne  system,  the 
rate  selector  switch  was  altered  to  allow  the  operator 
to  select  azimuth  rates  of  either  1,  5,  or  10  seconds 
per  revolution,  at  his  discretion. 

C.  The  elevation  drives  were  changed  to  yield  a  con¬ 
stant  step  of  5  degrees  in  elevation  at  the  comple¬ 
tion  of  each  azimuth  sweep.  In  the  automatic 
mode,  the  elevation  angle  for  the  first  revolution 
was  set  at  2-1/2  degrees  to  reduce  horizon  clutter 
induced  in  the  5-degree  field  of  view. 

The  resultant  sweep  pattern,  used  by  all  automatic  2ir 
scanners  during  the  OPAQUE  deployment  series,  is  illus¬ 
trated  in  Table  2.1.  In  this  pattern,  the  eighteen  azimuth 
revolutions  are  completed  in  180  sec  and  the  reset  and 
stop  intervals  are  automatically  deleted  from  the  recorded 
data. 


Table  2.1.  Automatic  2ir  Scanner  Sweep  ^tern 


Azimuth 
Revolution  No. 

Comtant 

Elevation 

An^ 

Upper  Hemiaphere 
Zenith  Angle 

Lower  ffcmiaphere 
Zenith  Angle 

! 

2.5* 

87.5* 

92.5* 

2 

7.5* 

82.5* 

97.5* 

3 

12.5* 

77.5* 

io2.r 

4 

17.5* 

72.5* 

107.5* 

S 

22.5* 

67.5* 

112.5* 

6 

27.5* 

62.5* 

117.5* 

7 

32.5* 

57.5* 

I22.r 

8 

37.5* 

52.5* 

127  5* 

9 

425* 

47.5* 

1315* 

10 

47.5* 

42.5* 

137.5* 

11 

1  52.5* 

37.5* 

142.5* 

12 

1  57.5* 

32.r 

147.5* 

13 

'  62.5* 

27.r 

152.5* 

U 

1  67.5* 

22.5* 

157.5* 

15 

1  72.5* 

17.5* 

162.5* 

16 

77.5* 

I2.r 

I67.r 

17 

82.5* 

75* 

l7Zr 

18 

87,5* 

2.5* 

177.5* 

3 


Gyro  Controi  Panel.  All  azimuth  and  elevation  ana¬ 
logs  generated  by  the  2ir  scanner  orientation  potentiome¬ 
ters  are  related  to  the  aircraft  frame.  If  the  aircraft  is 
Qying  straight  and  level,  all  scanner  orientation  data  is 
relatively  true.  Aircraft  bearing  is  recorded  on  the  project 
dato  logger  so  that  corrections  can  be  made  during  dau 
processing  for  bearing  variations  during  fli^t.  However, 
since  no  corrections  are  made  for  orientation  errors 
induced  by  pitch  or  roll,  a  reliable  monitoring  of  these 
maneuvers,  if  any,  was  desirable.  As  a  result,  a  Sperry 
system  SP-40  vertical  gyroscope  was  rigidly  mounted  to 
the  floor  of  the  aircraft  directly  below  the  upper  hem- 
sphere  27r  scanner.  Using  this  gyro  as  a  vertical  rderence, 
aircraft  pitch  and  roll  analogs  were  generated  and  fed  to 
the  in-flight  data  logger.  The  erection  and  control  circuits 
required  to  drive  and  monitor  the  project  gyro  are  con¬ 
tained  in  the  gyro  control  panel. 

The  optical  components  within  the  scanner  assembly 
are  illustrated  in  Fig.  2-5  which  has  been  compressed 
slightly  to  simplify  the  illustration.  The  flux  beam  emerg¬ 
ing  from  the  telescope's  exit  pupil  is  slightly  divergent 
such  that  an  image  approximately  IQmm  in  diameter  is 
incident  upon  the  detector’s  25mm  diameter  photo¬ 
cathode.  The  detector  assembly  which  has  been  described 
earlier  in  Duntley,  etai  (1975),  is  illustrated  in  Figs.  2-6 
and  2-7.  The  essential  characteristics  of  this  detector  cir¬ 
cuit  are  listed  below. 

Power  Input:  28VDC  at  approx.  300  ma  max. 

Signal  Output:  Selectable,  0  to  ±10VDC 

Dynamic  Range:  Calibrated  over  5  decade  change  in 
radiance 

The  output  of  the  scanner's  radiometer  circuit  was 
sampled  sixty  times  per  second  by  the  on-board  system 
data  logger.  For  most  applications,  these  constant  zenith 
angle  but  azimuthally  increasing  data  points  were  averaged 
to  yield  one  radiance  value  for  each  six  degrees  change  in 
azimuth.  Thus  the  sky  or  terrain  radiances  measured  by 
this  system  normally  ended  up  in  a  60x18  data  array 
wherein  the  2ir  steradians  were  represented  by  an  average 
radiance  value  every  six  degrees  in  azimuth  and  every  five 
degrees  in  elevation  iLe.  zenith  angle).  The  complete 
(unaveraged)  data  base  was  retained  in  back-up  for  diag¬ 
nostic  troubleshooting  and  special  purpose  displays. 


Fig.  2-i.  V»bility  Laboratory  integrated  radiation  detection 
asaembty. 


The  spectral  characteristics  of  the  scanner  system's 
radiometer  are  illustrated  in  Fig.  2-8  and  Table  2.2.  The 
detector,  an  EMR541E  series  photomultiplier  (S-20  spec¬ 
tral  response)  was  fitted  with  a  selection  of  interference 
and  absorption  filters  to  yield  the  illustrated  spectral 
responses.  The  relative  response  labeled  number  4  is 
similar  to  the  photopic  response  of  the  human  eye,  and 
thus  measurements  made  using  this  spectral  response 
represent  good  approximations  of  human  performance. 

Some  typical  atmospheric  properties  characterizing 
the  passbands  illustrated  in  Fig.  2-8  are  summarized  in 
Table  2.2.  The  last  row  in  Table  2.2  tabulates  the 
appropriate  properties  for  the  true  isotopic  response  for 
comparison  with  the  pseudo-photopic  properties  associated 
with  response  number  4.  It  is  clear  that  even  though 
responses  four  and  nine  have  substantially  different 
effeaive  passbands,  where  the  effective  passband  is  the 
area  under  the  normalized  rdative  spectral  response  curve 
and  is  equal  to  the  width  of  a  rectangular  passband  of 
equivalent  area,  the  atmospheric  properties  characterizing 
these  two  responses  are  quite  similar. 


(I)  elevation  prism 


(2)  AZIMUTH  PRISM 


(3)  ACHROMAT 


Fig.  2-5.  Automatic  2n  scanner  optica]  schematic.  (1)  Elevation  prism  #396732,  ieager,  41mm  x  lOimm  x  30mm.  (2)  Azimuth  prism 
#40995,  Edmunds,  36mm  x  25mm  x  25mm.  (3)  Achromat  doublet,  EFL*82,  B^*75.10,  FPL *83.25,  d|>*8.89.  (4)  Double-convex, 

EFL-53.64,  BFL-S2.04,  FFL-52.00. 


Fit*  2*7.  Visibility  Laboratory  radiation  detection  assembly  electrical  schematic. 


WAVELENGTH,  a  (nm) 

Of.  2-8.  Standard  spectral  responses  •  Pro^  OPAQUE. 


TeMc  2.2.  Sptctnl  ChwacierMia  Summiry  for  fnim  Of  AQUE 


3.  SKY  A  TEKKAIN  RADIANCE  DISPLAYS 

Partial  summaries  of  data  availabilities,  applications, 
and  selected  display  formats  related  to  the  airbme  meas* 
urements  condwned  with  the  C-130  mounted  system  have 
been  presented  previously  in  each  of  the  technical  reports 
listed  in  Appendix  C.  llie  most  comprehensive  of  these 
earlier  reports  is  presented  by  Duntiey,  etai  (1978c). 
For  the  nugority  oi  cases,  however,  these  summaries 
addressed  data  other  than  the  specific  radiance  arrays  pro¬ 
duced  by  the  Automatic  2v  Sc^er  Systems.  Therefore, 
in  the  following  sections  of  this  current  report  we  will 
address  specifically  the  sky  and  terrain  radiance  measure¬ 
ments  collected  during  a  series  of  European  flights  am- 
ducted  in  coigunction  with  program  OPAQUE,  Perm 
(1978)  and  Johnson,  etai  (1979). 

As  noted  in  Duntiey,  etai  (1978c),  there  were 
nearly  ninety  separate  flights  conducted  during  the  five 
different  OPAQUE  deployments.  Many  of  these  flights 
were  made  under  poor  weather  conditions  which  resulted 
in  only  partial  completion  of  the  anticipated  flight  profile. 
Thus,  not  every  fU^t  would  yield  sky  and  terrain  radiance 
measurements,  even  though  there  were  measurements  of 
volume  scattering  coefficient  profiles.  A  summary  of  the 
number  of  flights  yielding  each  of  the  two  basic  classes  of 
dau  i.e.  1)  sky  and  terrain  radiances,  and  2)  scattering 
coefficient  profiles,  is  shown  in  Table  3.1. 


TaMe  3.1.  Sumfinry  of  OPAQUE  Diu  Attempts 


Number  oT  Pfipits  | 

Depioynwnt 

Sky  a  Terrain 
Radiance 

Scattering  CoefScieni 
Profiles  Only 

Aborts 

Totals 

OPAQUE  1 

9 

0 

4 

13 

OPAQUE  11 

9 

4 

0 

13 

OPAQUE  III 

» 

3 

I 

13 

OPAQUE  IV 

14 

13 

2 

29 

OPAQUE  V 

16 

3 

1 

20 

Totals 

57 

23 

t 

U 

As  of  September  1981,  there  were  42  successful  and 
5  partially  successful  recoveries  of  sky  and  terrain  radiance 
sets,  where  a  "successful  recovery"  is  defined  as  a  flight 
yiel^ng  a  full  array  in  at  least  two  spectral  bands  at  two  or 
more  different  altitudes,  a  single  (2+2)  profile  in  the 
notation  of  Duntiey,  etai  (197te).  Lists  of  these 
recovered  data  are  tabulated  in  several  convenient  formats 
in  Section  4.  The  ten  flights  for  which  sky  and  terrain 
radiances  were  measured,  but  for  which  no  data  were 
recovered  are  documented  on  data  tapes  containing  severe 
parity  error  problems  which  preclude  their  immediate  pro¬ 
cessing.  The  displays  presented  in  the  remainder  of  this 
report  are  typical  of  those  available  from  the  successfully 
recovered  data  arrays. 

The  optimum  set  of  sky  and  terrain  measurements 
for  a  single  flight  results  in  a  pair  of  data  arrays  represent¬ 
ing  radiance  distributions  as  seen  from  each  of  four 
different  altitudes  and  in  each  of  four  different  spectral 
bands,  a  "Dual  (2+4)”  profile  in  earlier  notation.  Thus 
for  e^  optimum  flight  episode,  sixteen  pairs  of  data 


arrays  will  emerge  from  the  post  flight  data  processiiig 
sequences.  In  order  to  cope  with  the  analysis  and  evalua¬ 
tion  of  these  data,  several  different  graphical  formats  have 
been  developed,  which  are  illustrated  in  Figi.  3-2,  3-3  and 
3-4.  Eadt  cff  these  displays  may  be  derived  from  the  same 
basic  tabular  array  ehtch  contains  a  radiance  value  at 
every  six  degrees  in  azimuth  and  five  degrees  in  zenith 
angle.  It  is  possible  however,  to  retrieve  the  data  at 
snudler  azimuth  intervab  with  special  handling. 

As  an  aid  to  the  interpretation  of  radiance  distribu¬ 
tions  illustrated  in  Figs.  3-2,  3-3  and  3-4,  a  pictorial 
representation  of  the  measurement  field  is  shown  in 
Fig.  3-1.  This  artist's  conception  of  the  4ir  radiance  field 
which  surrounded  the  aircraft  is  annotated  with  the  angu¬ 
lar  notations  and  terminologies  used  to  define  the  data 
displays  which  follow. 


ZENTTH  angle  •  0* 


ZENITH  angle  •  l«r 


ng.  3-1.  Sky  and  terrain  coordinate  system. 

3.1  Four  Quadrant  •  Constant  Azimuth  Plots 

The  Four  (Quadrant  -  Constant  Azimuth  plots  illus¬ 
trated  in  Fig.  3-2  are  the  basic  diagnostic  displays  used  in 
the  evaluation  of  the  Automatic  2ir  Scaimer  ^u.  In  each 
plot  of  this  type,  the  azimuth  coordinates  have  been 
adjusted  to  index  relative  to  the  sun  position  at  the  time 
of  the  measurement.  Thus,  the  data  points  coded  as 
azimuth  zero  all  were  measured  as  the  Smner’s  routing 
field  of  view  intercepted  the  vertical  sun-zenith  plane, 
looking  toward  the  sun.  Convmely,  the  dau  points  coded 
as  azimuth  180*,  were  all  measured  as  the  field  of  view 
intercepted  the  same  plane  looking  directly  away  from  the 
sun.  The  two  sets  of  points  coded  azimuth  90*  and  270* 
similarly  represent  the  measurements  in  the  two  cross-sun 
directions.  In  all  cases,  a  zenith  angle  of  zero  specifies  the 
point  in  the  sky  directly  overhead,  a  zenith  aitgle  of  90* 
represents  the  horizon,  and  a  zenith  angle  of  180* 
represents  the  point  on  the  terrain  directly  below. 
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60  80  100  120  140  160  180 

ZENITH  ANCLE  (Defiees) 

FIf.  3-2.  Four  Quadrant  •  constant  azimuth  radiance  scans  from 
filter  2  of  Flight 

The  two  sets  of  curves  in  Fig.  3-2  illustrate  several 
nuuor  atmospheric  characteristics  one  consistently 
observes  in  analyzing  data  of  this  iort,  and  in  relating  the 
measurements  to  the  visually  observed  world. 

a.  The  sky  nuliances  are  orderly  in  the  magnitudes. 
The  up-sun,  azimuth  zero  data  are  the  highest,  the 
down-sun  azimuth  180"  the  lowest,  and  the  cross- 
sun,  azimuths  90"  and  270"  fall  neatly  in  between. 

b.  The  sky  radiances  away  from  the  sun  are  substan¬ 
tially  higher  at  the  lower  altitude  (264m  AGL)  as 
one  would  expect  from  the  increased  path  radiance 
influence. 

c.  The  horizon  radiances  are  substantially  higher  than 
either  the  zenith  or  nadir  radiance  at  the  higher 
Ol03m  AOL)  altitude. 

d.  The  variations  in  terrain  radiance  are  smoothed 
markedly  as  the  altitude  of  measurement  changes 
from  264  to  3103m  AGL. 

e.  The  high  radiances  near  the  sun  drive  the  radiome¬ 
ter  and  the  measurements  show  evidences 

of  being  adversely  influenced  by  stray  light  within 
the  optical  system. 
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S.2  Site  AteMtttar  PMt 

A  special  set  of  measurements  were  normally 
included  in  each  Automatic  2ir  Scanner  sequence  in  an 


AZIMUTH  ANGLE  (Degrees) 

Fig.  3-3.  Solar  almucantar  scans,  direct  and  composite.  Flight 
C-468,  Alter  2,  altitude  --  2960m  AGL. 


attempt  to  more  fully  document  the  radiances  in  the  near 
vicinity  of  the  sun.  ITiis  measurement  sequence  referred 
to  as  the  SUN  MODE  sequence,  involved  a  rescan  of 
three  selected  azimuth  revolutions  after  first  inserting  a 
3  log  neutral  density  filter  into  the  scanners  optical  path. 
The  insertion  of  this  filter  prevented  the  near-sun  radi¬ 
ances  from  saturating  the  radiometer  circuit  and  permitted 
the  measurements  illustrated  in  the  top  graph  of  f^.  3-3. 

The  measurements  made  through  the  neutral  den¬ 
sity  filter  were  referred  to  as  the  SUN  MODE  REVOLU¬ 
TIONS  and  their  plot  (3-3t)  is  annotated  SUN  REV  for 
identification.  One  notes  from  the  plot  that  all  data  points 


ZENITH  ANGLE  (Degrees) 


ZENITH  ANGLE  (D^rees) 


Fig.  3-4.  Composite  sky-terrain  radiance  scans.  Flight  C-468. 
filter  2. 

lie  above  the  corrected  minimum  calibratable  value 
(CMCV)  measurable  by  the  5  decade  radiometer  span. 
The  "corrected  minimum  calibratable  value"  (CMCV) 
represents  that  radiance  which  is  the  minimum  radiance 
value  that  the  detector  can  usefully  discriminate  after 
being  properly  corrected  for  spectral  band  and  sensitivity 
shifts.  For  these  systems  it  is  the  radiance  level  at  which 
the  detector  output  is  either  three  decades  below  the  radi¬ 
ance  established  by  the  systems  internal  Isolite  luminous 
reference  cell  or  at  which  the  signal  to  noise  ratio  becomes 
unacceptable,  whichever  occurs  first. 

Measurements  made  without  the  neutral  density 
filter,  made  immediately  preceding  the  SUN  MODE,  and 
for  the  same  three  scanner  revolutions  are  illustrated  in 
Fig.  3-3b  and  are  identified  as  SKY  REV  data.  These 
measurements,  as  illustrated  in  both  Figs.  3-3b  and  3-2 
exceed  the  corrected  zero  scale  value  (CZSV)  of  the 
radiometer  in  the  region  near  the  sun,  and  thus  are  not 
reliably  retrievable.  The  "corrected  zero  scale  value" 
(CZSV)  represents  that  radiance  which  is  the  maximum 
radiance  value  that  the  detector  can  usefully  discriminate 
after  being  properly  coirected  for  spectral  band  and  sensi¬ 


tivity  sliifts.  F6r  these  systems  it  Is  optically  a^usted  to 
be  ^  radiaiice  level  which  oocuis  two  decMles  above  the 
radiance  esublished  by  the  sysiem*s  internal  bolite  lumi¬ 
nous  reference  cell.  The  reference  radianoe  is  establtihed 
such  that  the  system’s  five  decade  span  will  cover  aU  radi¬ 
ance  levels  nomudly  encountered  during  a  field  measure¬ 
ment  program. 

A  computer  generated  composite  of  the  two  data 
sets  discussed  above,  Le,  the  SUN  REV  and  SKY  REV 
data,  is  illustrated  in  Fig.  3-3c  and  is  identified  as  the 
EDIT  REV  dau.  In  this  composite  set,  the  computer  has 
matched  those  dau  points  where  the  SUN  REV  and  SKY 
REV  values  were  equal  after  adiusting  for  filter  ofbets, 
and  combined  this  matched  set  into  a  single  array  in  which 
the  values  from  the  SUN  REV  are  used  in  those  regions 
where  the  SKY  REV  data  were  beyond  the  CZSV,  i.e. 
exceeded  the  calibrated  span  of  the  radiometer.  Thus,  the 
inclusion  of  these  EDIT  REV  dau  points  in  the  final  daU 
array  representing  the  measured  values  sky  radiance 
results  in  a  single  array  containing  a  measured  value  for 
every  point  in  the  sky,  including  the  immediate  vicinity  of 
the  sun. 

As  with  the  interpreution  of  Fig.  3-2,  there  are 
several  repeating  characteristics  one  regularly  observes 
when  analyzing  the  displays  of  Fig.  3-3. 

a.  The  average  radianoe  over  the  scanner’s  five  degree 
circular  field  is  sudi  that  the  peak  values  in  the  SUN 
REV  are  normally  about  3.5  decades  above  the  aver¬ 
age  sky  radiances. 

b.  Since  the  SUN  REV  displays  plot  dau  representing 
the  fixed  azimuth  scan  vi^i^  occurs  nearest  the  true 
solar  zenith  angle  plus  dau  from  the  preceding  and 
following  scans,  one  nuy  deduce  the  extent  to 
which  the  solar  disc  was  centered  in  the  field  of  view 
during  the  nominal  solar  intercept.  The  key  to  this 
deduction  is  the  symmetry  of  the  before  and  after 
scans.  That  is,  if  the  preening  and  following  almu- 
canur  scans  are  exactly  congruent,  then  they  were 
symmetrically  located  above  and  below  the  disc  and 
the  central  scan  would  fully  capture  the  disc  within 
its  field  of  view.  There  is  of  course,  due  to  the  S" 
field  of  view  of  the  instrument,  a  family  of  increas¬ 
ing  assymetries  for  which  the  disc  is  still  captured  by 
the  SUN  REV  but  with  decreasing  certainty. 

c.  Comparison  of  these  almucantar  radiances  with 
modelled  aureole  and  surround  radiances  indicate 
the  measured  near-sun  sky  radiances  may  be  too 
high. 

3.3  Compaslte  Sky-Terrain 
Radianoe  Scans 

The  sky  and  terrain  radianoe  displays,  illustrated  in 
Fig.  3-4  are  a  condensed  version  of  the  I^r  Quadrant 
plot  shown  in  the  preceding  Fig.  3-2,  and  are  created  from 
the  same  ubular  dau  base.  In  these  composite  plots,  the 
radiance  variations  throughout  the  sun-zenith  plane  in 
both  the  up-sun  and  down-sun  directions  are  plotted  as  a 
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single  curve,  and  labeled  *PH[-0  4  •IKT.  Similarly,  tbe 
dau  in  the  croas-aun  plane  are  plotted  oondnuoualy  and 
labeled  "PHI  *  90  4  -270”.  In  all  cases,  a  aenith  view 
angle  of  0*  represents  the  point  in  the  si^  diiectly  over¬ 
head,  and  a  zenith  view  angle  of  ISO*  represents  the  point 
on  the  underlying  terrain  directly  below  as  illustrated  in 
Fig.  3-1.  The  negative  signs  associated  with  angular 
indices  are  an  artifact  of  the  computer  sorting  tedmique, 
and  should  not  be  rigorously  interpreted.  In  general,  the 
negative  sign  is  used  merely  to  dedgnate  the  two  different 
azimuthal  directions  associated  with  each  spherical  scan. 

The  utility  of  the  composite  plots  shown  in  Fig.  3-4 
lies  in  the  compact  visualization  of  the  radiance  field  sur¬ 
rounding  the  aircraft  at  the  time  of  measurement.  The 
radiance  field  characteristics  discussed  briefly  in  conjunc¬ 
tion  with  the  Four  Quadrant  displays  tend  to  be  more 
clearly  and  succinctly  defined  in  the  composite  plots. 
Thus,  visual  comparisons  by  the  analyst  tend  to  be  simpler 
and  faster.  Having  the  entire  upper  hemisphere  radiance 
pattern  presented  symmetrically  about  the  center  of  the 
plot  is  intuitively  satisfying  to  most  users  and  enhances 
ones  association  of  the  data  display  with  ones  concept  of 
the  real  world. 


4.  DATA  IVALUATKm 

The  data  displays  ittustrated  in  Pig.  3-4  repraasnt  the 
fonnat  used  on  all  of  the  Automatic  2w  Scanner  meeroie- 
ments  that  were  retrievahle  from  the  (VAQUB  aeiiae  of 
deployments.  Thus  a  plot  siniilar  to  3-4a  is  available  for 
each  of  a  broad  aelectkm  of  measurement  acenarioa.  The 
situatioos  for  which  these  aky  and  terrain  radiance 
representations  have  been  proceseed  are  aummaiized  in 
Table  4.1,  segregated  by  geograiducal  location  in  Table  4.2, 
segregated  by  "number  of  ploCs  per  flight"  in  Table  4.3, 
and  segregate  by  upper  hemisphm  doud  dassificatioQ  in 
Table  4.4.  As  indimted  in  these  summarizing  tablet,  the 
prepondefince  of  the  radiance  data  were  collected  over 
five  of  the  OPAQUE  stations:  Meppen  and  Biikhof,  Ger¬ 
many,  Rodhy,  Denmark,  Yeovilton,  England  and  Bruz, 
France.  It  is  also  apparent  that  there  is  a  good  distributkm 
of  reiattvely  complete  sets  /.e.  all  four  spectral  bands  at  all 
four  altitudes  (9  flights)  or  all  four  spectral  bands  at  each 
of  two  altitudes  (15  flights).  Thus,  with  a  data  base  of 
nearly  500  graphical  displays,  one  may  confidently  address 
the  use  of  these  date  in  a  varidy  of  technical  applications, 
dependent  only  upon  an  assessment  of  the  reliability  and 
qudity  of  the  measurements. 


TaMe  4.1.  Summary  of  Sky  and  Terraifi  Sadianoe  Meaauramenu 
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N«a;  Hie  neminel  ttgit  Utitudei  Heied  in  cdumni  4  tlirougi  7.  the  akitudei  amorieted  with  each  of  the  Iher  2  meaemementa.  ere  coded  in  cotmm  t  dwomh  II 
to  deeltnaie  thoie  ahitudee  at  which  additional  ipactm  data  are  available.  Thus  for  flgH  371,  the  entriei  *Aa~”  in  oohimm  I  throiigi  U  Mteam  thot  meooum- 
fflema  in  each  of  the  four  ipecttal  band  ware  made  ai  each  of  two  fight  ahitudee.  A*2t2ni  and  S«17Siin. 
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TflUc  4.4.  RetrievMl  py^  Sorted  by  Upper  llemiapliefe  Cloud  niwiflrrtinn 
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Note  1:  The  cfoud  dassiflcition  dtte  listed  in  T«Me  4.4  have  been  extracted  from  ftve  previous  reports  rdaftcd  to  the  OPAQUE  fll|ht  proffam,  Dunilep,  0  oL  (19T7. 
I97ta,  197th)  and  Johnson  and  Gordon  <1979  and  19M). 

Note  2:  Since  many  of  these  lli|his  were  nude  during  rapidly  charing  weather  condiitons,  substantial  subjectivity  is  involved  in  the  assignment  of  any  given  fliglit  to  a 
specific  doud  condition. 

Note  3:  Flight  type  notation  i.e.  D(24-4)  as  indicated  in  previous  reports  b  interpreted  as  a  Dual  profile  each  containing  roeasuiemenu  in  two  spectral  bands  at  each  of 
four  diffi^i  altitudes. 


4.1  Stray  Light  ErMciim 

An  evaluation  of  the  reliability  of  the  sky  and  ter¬ 
rain  radiance  measurements  obtained  by  the  Autontatic  2ir 
Scanner  systems  has  been  conducted  along  several 
different  avenues  throughout  the  measurement  program. 
In  general,  the  directional  orientation  and  radiometric  cali¬ 
bration  relationships  have  been  found  to  be  reproducible 
and  within  design  specifications,  /.e.  ±0.5*"  in  azimuth 
and  elevation,  and  ±3%  in  absolute  radiance.  However, 
as  one  might  suspect  for  an  optical  scanner  of  this  design, 
see  Fig.  2-5,  there  were  evidences  in  the  dato  indicating 
significant  stray  light  influences  as  the  systems  field  of 
view  approacheid  the  location  of  an  unobacured  solar  disc. 

The  earliest  hints  of  potential  uncertainty  in  the  sky 
radiance  dau  appeared  in  comparisons  between  values  of 
integrated  upper  hemisphere  scanner  radiances  and  the 
dual  irradiometer  measurements  of  irradiance.  These 
comparisons  of  two  independently  determined  values  of 
irradiance  illustrated  an  irregular  and  seemingly  non- 
systematic  mismatch  in  the  two  irradiance  d^erminations. 
Several  possible  explanations  related  to  sampling  intervals. 


optical  averaging,  and  irregular  cloud  conditions  were  dis¬ 
cussed  at  length,  but  none  seemed  adequate  to  identify 
the  variances  In  the  data.  As  a  result  of  this  uncertainty, 
computations  of  atmospheric  optical  properties,  such  as 
those  reported  in  Duntley  era/.  (1976),  used  irradtances 
derived  from  scanner  arrays  as  the  primary  data  input 
since  they  were  normally  used  in  ratio  form  with  their 
companion  radiance  data,  and  thus  instrumental  uncertain¬ 
ties  would  in  most  cases  be  minimized,  liowevw,  there 
was  a  continuing  concern  that  there  was  an  undetected 
error  somewhere  in  the  experiments  procedure. 

With  the  advent  of  the  ffering  nnodel  for  estimating 
atmospheric  contrast  transmittance,  Hering  (1981),  an 
operationally  efficient  mechanism  was  piDVided  for  diag¬ 
nostic  analysis  of  atmospheric  radiance  distributions. 
Comparisons  of  sky  and  terrain  redianoe  distributioos  cal¬ 
culated  by  the  model  to  those  measured  by  the  airttome 
scanner  systetiw  are  consistently  excellent  in  all  regions 
except  near  the  sun.  Figures  4-1  and  4-2  which  have  been 
reproduced  from  Hering  (1981)  are  typical  of  the  model- 
measuremMit  comparison.  It  must  be  noted  however  that 
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Fig.  4-1.  Comperison  of  measured  and  calculated  sky  and  terrain 
radiance  distributions  at  200m  (Hering,  1901). 


Fig.  4-2.  Comparison  of  measured  and  calculated  sky  and  terrain 
radiance  distributions  at  6000m  (Hering.  1981). 


the  measured  data  in  Figs.  4-1  and  4-2  have  been  deleted 
in  the  region  near  the  solar  zenith  angle  of  43’’.  The  need 
for  these  deletions,  and  the  source  of  the  earlier  irradianoe 
inconsistencies  were  both  immediately  illustrated  by  the 
model  comparison.  When  model  calculations  such  as 
those  shown  in  Fig.  4-2  were  compared  with  sky  radiance 
measurements  as  illustrated  in  Fg.  3-4,  the  near  sun 
measurements  were  consistently  too  high  with  respect  to 
the  model  calculation.  Since  modelled  and  measured 
values  were  consistently  equivalent  in  all  other  portions  of 
the  scene,  the  near  sun  radiance  measurements  were  con¬ 
sidered  to  be  the  primary  source  of  the  previously 
unidentified  experimenul  error. 


4.2  ExparlMital  ValMatlaa 

With  both  irradianoe  measurements  and  model 
predictions  indicatini  erroneous  measurements  of  hiear 
sun  radiances,  an  experimental  measurement  sequence 
was  devised  as  a  test.  The  procedure  was  straightforward. 
The  Automatic  2ir  Scanner  SyMm  was  recalibrated  and 
set  up  for  operation  with  manual  tracking  control  and  with 


the  radiometen  sensitivity  adjusted  to  permit  direct  solar 
viewing.  The  following  measurement  sequence  was  then 
repeated  several  times  daily  for  a  period  of  about  ten  days. 

A.  Measure  apparent  radiance  of  sky  at  nine  positions 

along  the  solar  almucantar  starting  at  ^  and 
ending  at  ^  +  55  ”,  where  ^  =  azimuth  angle 

of  the  center  of  the  scanner’s  field  of  view  and 

=  azimuth  angle  of  the  sun. 

B.  Repeat  eadi  measurement  in  step  A  with  the  solar 
disc  obscured  by  a  4”  opaque  mask. 

C.  Compare  the  radiance  measurements  made  at  each 
angular  position  with  and  without  the  occulting 
mask  to  determine  the  off  axis  stray  light  com¬ 
ponent  in  each  pair  of  measurements. 

The  experimental  procedure  outlined  above  verified 
the  existence  of  a  substantial  stray  li^t  contribution  to  the 
measured  signal  at  all  angles  within  30”  from  the  sun. 

The  dau  were  collected  at  the  surface  (-'300ft 
MSL)  on  days  with  varying  aerosol  loads  ranging  from 
heavy  haze  to  relatively  dear. 

As  mentioned  in  step  A  of  the  procedure  outlined 
above,  measurements  of  itoucantar  radiance  were  made 
at  nine  angular  positions  with  respect  to  the  solar  azimuth. 
They  began  at  =  0”  and  then  advanced  to  5”,  7”, 
10®,  15®,  20®,  25®,  30®,  and  55®.  Each  measurement 
sequence  took  approximately  ten  minutes  to  complete, 
during  which  there  were  2  to  3  degrees  change  in  solar 
azimuth.  Thus,  since  the  measurement  sequence 
advaimd  in  the  same  direclion  as  the  sun's  transit,  the 
measured  radiances  occurred  slightly  closer  to  the  solar 
disc  than  indicated  by  the  nominal  locations  indicated 
along  the  horizontal  axis  of  Fig.  4-3.  The  measurements 
made  at  0®  are  not  illustrated  since  the  hardware  associ¬ 
ated  with  the  occulting  procedure  was  inadequate  to  repro- 
ducibly  locate  a  small  disc  accurately  at  the  center  of  the 
scanner's  field  of  view  and  thus  only  the  unooculted  meas¬ 
urements  were  obtained.  One  does  note  however,  in 
colunm  5  of  Table  4.5,  the  progressive  decrease  in  the 
average  field  radiance  at  ^^,-0*as  the  aerosol  loading 
gradually  increases.  The  55®  dau  are  not  illustrated  since 
at  this  angular  distance  from  the  sun  there  were  never  any 
evidences  of  stray  tight,  and  the  indusson  of  the  dau 
points  would  needlessly  expand  the  scale  of  the  illustra¬ 
tion. 


Dau  representing  the  unooculted  to  occulted  radi- 
anoe  ratios  at  each  of  the  remainmg  seven  almucantar 
positions  for  each  of  ten  separate  almucantar  scans  are 
plotted  in  Fig.  4-3.  In  this  figure,  each  of  the  ten  plots 
represents  the  magnitude  of  the  scanner's  stray  tight  error 
for  a  different  aerosol  condition  as  noted  in  Table  4.5. 
The  stray  light  ratios,  A$l,  which  are  indicated  along  the 
vertical  axis  of  Fig.  4-3  are  defined  by 
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Flt.4>3.  Experimontally  determinod  straytight  effecu. 


where  for  the  denominator  measurement  the  scanner 
entrance  aperture  is  centered  in  the  shadow  cast  hy  an 
opaque  circular  disc  which  subtends  a  field  of  4*  as  seen 
from  the  scanner  entrance. 

It  is  apparent  from  the  dau  illustrated  in  Fig.  4-3 
that  the  off-axis  light  entering  the  scanner’s  optical  system 
is  substantial  at  all  viewing  angles  near  the  sun,  and 
increases  markedly  both  as  the  fidd  of  view  approadies 
the  solar  disc  and  as  the  atmospheric  aerosol  loading 
decreases. 

The  influence  of  the  aerosol  load  on  stray  light  error 
and  thus  the  Adefity  of  the  scanner’s  radiance  measure¬ 
ments  is  illustrated  best  perhaps  by  considering  two  con¬ 
trasting  examples.  In  both  examples  consider  the  meas¬ 
urement  of  sky  radianoe  at  the  same  distance  from  the 


sun,  however  in  the  first  example  the  measurement  is 
made  with  the  scanner  at  6  km  in  altitude  and  in  the 
second  at  only  1  km. 

At  6  km  the  aerosol  loading  is  typically  very  low 
with  atmospheric  optical  scattering  ratios  Le.  total 
scattering/Rayleigh  scattering,  Hering  (1981),  in  the  order 
of  1.2  to  1.5  whence  the  apparent  radianoe  of  the  predom¬ 
inantly  blue  sky  is  relatively  low.  Conversely,  the  direct 
solar  irradiance  falling  obliquely  upon  the  scanner’s 
entrance  aperture  is  only  slightly  attenuated  by  its  passage 
through  the  thin  upper  atmosphere  and  thus  gives  rise  to 
relatively  higli  level  off-axis  reflections  within  the  optical 
system.  These  conditions  would  give  rise  to  high  stray 
light  ratios,  such  as  those  Illustrated  by  the  curves 
from  tests  2,  6  and  10  shown  in  Fig.  4-3. 

In  the  example  with  the  scanner  at  only  1  km  alti¬ 
tude  the  measurement  of  apparent  sky  radial  is  made 
through  a  mudi  more  turbid  medium  in  which  the  typi¬ 
cally  heavier  aerosol  load  generates  substantially  hi^ier 
optical  scattering  ratios  and  their  resultant  higher  propor¬ 
tion  of  forward  scattering.  In  this  example  the  apparmt 
radianoe  of  the  now  hazy  sky  is  high  relative  to  that  meas¬ 
ured  at  6  km,  however  the  same  scattering  mechanism 
a4iich  has  increased  the  apparent  sky  radianoe,  has  also 
contributed  to  the  attenuation  of  the  direct  solar  irradiance 
falling  upon  the  scanner’s  aperture  and  thus  reduced  the 
intensity  of  the  system’s  internal  stray  light.  This  second 
example  illustrates  the  mit^ating  influent  of  the  aerosol 
loading  upon  the  scanner  measurement’s  stray  fight  com¬ 
ponent.  These  low  altitude  conditions  would  give  rise  to 
lower  stray  light  ratios  such  as  those  illustrated  by  the 
curves  from  test  numbers  8,  1,  and  4  ^iown  In  Fig.  4-3. 

In  the  full  overcast  condition,  where  the  trradianoe 
from  a  sky  segment  just  off-axis  is  not  substantially  higher 
than  that  from  an  on-axis  segment,  or  at  viewing  angles  in 
excess  of  3fl*  from  the  solar  disc  there  are  no  discmible 
stray  light  influences  uton  the  scanner’s  radianoe  meas¬ 
urements  and  the  values  of  approach  one  as  illustrated 
along  the  lower,  right  hand  boundary  of  Fig.  4-3. 
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The  impUciaioiis  to  be  derived  from  these  data  are 
that  high  attitude  or  dear  day  measurementa  of  sky  radi- 
anoes  within  10*  of  the  solar  disc  win  have  stray  light 
ratios  of  as  much  as  20  to  35,  whereas  low  altitude  near¬ 
sun  measurements  through  more  turbid  atmospheres  will 
have  substantially  lower  stray  light  latios  of  about  3  to  6  at 
the  same  angular  distance  from  the  disc.  In  either  case  It 
is  dear  that  with  stray  light  infiuenoes  of  the  magnitude 
illustrated  in  Fig.  4-3,  the  dau  displays  such  as  Fig.  3-2a, 
in  order  to  be  fully  reliable,  should  in  fact  be  presented  in 
a  form  such  as  that  illustrated  in  Fig.  4-4.  (Jdhnaon  and 
Hering  19S1).  Fortunately,  the  same  modd  that  pin¬ 
pointed  the  erroneous  measurenmts  provides  the  vehide 
for  a  realistic  correction  procedure. 


na.4-4.  I^teasured  sky  and  terrain  radiance  distributions  for 
Flight  C-466,  filter  2.  2953m  AGL 


4.3  Stray  Light  Correctioits 

Inasmuch  as  the  data  in  Fig.  4-3  illustrate  that  the 
validity  of  near  sun  radiance  measurements  is  highly 
infhi^iced  by  the  variations  in  atmospheric  aerosol  load¬ 
ing,  the  seleaion  of  an  appropriate  corrective  procedure  is 
somewhat  limited.  Only  two  have  been  conddered  at  this 
point  and  they  are  described  briefly  in  the  following  para¬ 
graphs. 

The  first  procedure  is  to  use  an  extrapolation  from 
the  last  valid  measurement  nearest  the  sun  to  mlculate 
new  sky  radiance  values.  A  technique  used  by  Barteneva 
(1%0)  to  extrapolate  measured  directional  scattering  func¬ 
tions  from  16*  scattering  angle  to  0*  scattering  angle  was 
applied  to  the  ffight  C-466  sky  radiance  data  as  a  test.  In 
this  technique  one  assumes  the  log  of  the  sky  radiance  to 
be  a  linear  function  of  the  scattering  angle  cosine  between 
30*  and  CT  scattering  angle.  This  relatively  familiar  rela¬ 
tionship  is  expressed  as  an  assumed  equality  in  Eq.  (4.1), 

In  Ll  itAfi)  -  Cosfi  (4.1) 

where  is  the  sky  radtanoe  U.  upward  looking 

path  radianoe,  at  altitude  z  at  zenith  angle  s  and  at  scatter¬ 
ing  angle  B  from  the  sun.  The  technique,  is  reasonably 
effoctive,  even  though  it  tends  to  significuitly  underesti¬ 
mate  the  diffraction  peak  of  the  phase  function  as  illus¬ 


trated  in  Rg.  4-5a.  Fortunately,  when  uaed  to  tt  In  llea^ 
sun  sky  radiance  vahiea  wlucii  wOl  be  buegmted  in  eoo- 
junction  with  a  separme  apparent  sdar  disc  radiance  aa  die- 
cusaed  in  Duntlay,  arot  (1975)  the  residual  error  u  the 
total  IrridiaDoe  level  is  sm^.  For  example,  when  appBad 
to  die  fight  C-466,  47tem,  higli  altitude  data,  the  calctt- 
lated  irradianoes,  which  had  originally  been  aa  much  m 
14%  higher  than  thoae  maasuied  by  the  irtadioiiieier  sya- 
tem,  were  brought  to  within  ±4%  of  the  measured  values. 

A  second  technkpie  is  based  upon  the  demonstrsted 
ability  of  the  Hering  model  to  recreate  the  measured 
values  of  the  diffuse  radianoe  field  end  the  vertical  struc¬ 
ture  of  the  trradianoe  discrtbutfons.  In  this  procedure,  the 
valid  measufemeots  within  the  array  /.c.  die  radianoe 
measurements  more  than  30*  from  the  sun,  in  ooidunc- 
tion  with  the  moM  algorithms,  establish  both  the  form 
and  the  magnitude  of  corrected  radtanoe  Add.  The  pro¬ 
cedural  concept  assumes  that  the  single  scsttering  phase 
function  can  be  represented  analytically  with  good  preci¬ 
sion  by  two  term  Henyey-Greennein  functions  (I^ne, 
1968)  even  though  these  two  term  functions  also  tend  to 
underpredict  the  forward  10*  of  the  phase  function  as  is 
illustrated  in  Fig.  4-5b  which  has  been  adapted  from  Kat- 
tawar  (1975),  and  that  these  analytic  expressions  are 
directly  relatable  to  the  measured  optical  scattering  ratio. 
Johnson  et  al.  (1979)  and  Hering  (1981). 

Even  though  it  is  more  complex  computationally, 
there  are  specific  advantages  to  be  gained  in  pursuing  the 
further  development  of  this  second  technique.  Not  only 
do  the  Henyey-Greenstein  functions  more  dosely  approxi¬ 
mate  the  forward  scattering  regions  of  the  atmospheric 
phase  function  as  illustrated  in  Fig.  4-4b,  hut  the  ftexibili- 
ties  inherent  in  conditioning  the  Hering  model  inputs  to 
the  measured  environment  add  substantial  realism  to  the 
computed  radiances.  As  indicated  in  Appendix  B, 
abstracted  from  Hering  (1981)  and  which  summarizes  the 
operation  of  his  model,  the  nuxlel  inputs  are  designed  to 
take  full  advantage  of  the  most  descriptive  dau  available. 
The  key  inputs  are  the  optical  scattering  mixing  ratio, 
Qiz),  and  the  single  scattering  albedOv  «i(z),  where  Qiz)  is 
defined  as  the  ratio  of  the  toul  volume  scattering 
coefficient  at  altitude  z,  to  the  molecular  (or  Rayletgh) 
volume  scattering  coefficient  at  the  same  altitude  z;  and 
w(z)  is  defined  as  the  ratio  of  the  total  volume  scattering 
coefficient  at  altitude  z.  to  the  volume  attenuation 
coefficient  as  determined  at  altitude  z.  The  model  will 
from  these  ^lecified  inputs  compute  the  Henyey- 
Greenstein  functions  associated  with  the  input  Qiz)  values 
and  proceed  with  the  oompuutioo  of  sky  ridlaooe.  How¬ 
ever,  to  the  extent  that  more  complete  informetion  is 
available  to  define  Fiz^)  the  atmospheric  phase  fimctkm 
at  altitude  z,  the  model  will  reedity  accept  this  infonnatioii 
in  lieu  the  computed  values  esrimated  from  the 
specification  of  Qiz), 

The  fidelity  with  which  the  Hering  model  has  repro¬ 
duced  sky  radiance  meaturemeou  in  reglona  free  ftom 
stray  light  conummatkm,  «  illusirated  in  the  prelimincry 
run  shown  in  Figs.  4-1  and  4-2,  indicetee  the  ffigh  prabe- 
biltty  that  it  will  produce  accurate  eftimgtes  of  the  tnie  sky 
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Extnpolttion  assuming  scattering  phase  function  liriear  with  cosine  of 
scattering  angle. 


e 

Apprwumation  assuming  scattering  phaae  function  represented  by  two  term 
Henyey-Oreenstein  function  (from  Kaitawar,  (197S)]. 


Fig.  4-5.  Comptrison  of  soriterint  plMse  function  approxinuitfons 
with  modelled  cekulationg. 


radiance  in  most  of  those  regions  currently  contaminated. 
Thus  the  proposed  corrective  procedure  would  consist  of 
the  following  steps  for  each  set  of  radiances  one  might 
desire.  1)  From  the  experimental  measurements,  deter¬ 
mine  and  enter  as  input  the  vertical  profiles  of 


0(z),  and  a»U)  including  estimates  of  these  proHes 
in  the  relatively  dear  atmosphere  above  the  higliest  alti¬ 
tude  sampled  by  the  instrumented  aircraft  (nonmUy 
6  km),  2)  Let  the  model  compute  the  sky  radiances  in  the 
planes  illustrated  in  Fig.  3-4.  3)  If  necessary,  atfiust  the 
input  parameters  to  insure  an  acceptable  model  fit  to  the 
measurements  in  the  regions  between  30*  and  55*  from 
the  solar  disc.  4)  Substitute  the  values  of  sky  radiance  as 
calculated  by  the  model  into  the  measured  array  as  the 
approximated  correct  values. 

It  is  apparent  that  even  with  this  technique  there 
will  be  some  uncertainty  in  the  reliability  of  the  calculated 
values  of  sky  radiance  as  the  values  approach  the  disc, 
particularly  as  they  move  inside  IS*  where  the  Henyey- 
Greenstein  functions  begin  to  fall  short  of  the  Mie  predic¬ 
tions.  It  would  be  a  desireable  facet  o(  the  experimental 
validation  of  this  process  to  conduct  a  series  of  simultane¬ 
ous  measurements  using  the  automatic  2ir  scanner  as 
described  in  Section  2  and  the  4'  field  of  view  Contrast 
Reduction  Meter  (CRM)  described  in  Duntley,  etal. 
(1972).  These  CRM  data  would  provide  significantly 
more  reliable  near-sun  radiance  measurements,  particu¬ 
larly  in  the  difficult  solar  aureole  regions  Le.  at  scattering 
angles  between  zero  and  about  five  degrees,  against  which 
both  the  scanner  and  the  model  performances  could  be 
judged. 

5.  SELECTED  SAMPLE  DATA 

It  is  beyond  the  scope  of  this  report  to  include  the 
nearly  500  graphs  and  tables  which  represent  the  sky  and 
terrain  radiances  retrieved  from  the  OPAQUE  flight  series. 
It  is  desireable  however  to  illustrate  to  a  limited  degree  the 
variety  of  characteristics  inherent  within  the  dau  set. 
Therefore  the  figures  contained  in  this  section  should  be 
considered  as  an  hors  d'oeuvre  is  considered,  with  respect 
to  that  which  might  follow. 

Three  flights  have  been  selected  to  provide  this 
sample,  flight  C-392  from  the  OPAQUE  II  deployment 
during  the  Fall  of  1976  as  reported  in  Duntley  etai. 
(1978a),  flight  C-410  from  the  OPAQUE  III  deployment 
during  the  Summer  of  1977  as  reported  in  Duntley  ei  al 
(1978b),  and  flight  C-468  from  the  OPAQUE  V  deploy¬ 
ment  during  the  Summer  of  1978  as  reported  in  Johnson 
and  Gordon  (1980).  As  a  convenience  to  the  reader,  the 
flight  description  sheets  for  each  of  these  sample  flights 
have  been  extracted  from  the  parent  reports  and  are 
included  herein. 


THE  DATA  PRESENTED  IN  THE  FOLLOWING 
HGURES  HAVE  NOT  BEEN  CORRECTED  FOR 
THE  NEAR  SUN  STRAY  LIGHT  EFFECTS  DIS- 
CUSSED  IN  SECTION  4. _ 


Flight  C-392  is  one  of  the  fifteen  listed  in  Table  4.3 
as  having  eight  retrievable  data  arrays.  In  most  cases,  as 
illustrated  by  flight  C-392,  these  dght  arrays  represent 
measurements  made  in  eadt  of  four  ^)ectral  bands  at  two 
different  altitudes,  a  dual  (2+2)  profile  in  the  noution  of 
Duntley  era/.  (1978a).  The  feature  most  noticeable  in 
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this  series  of  displays  is  the  intermittent  appearance  of  the 
solar  peak.  In  order  to  easily  calculate  the  contrast 
transmittance  betiveen  the  two  altitudes  shown  in  Fig.  S-la 
and  S*lb  for  example,  the  radiance  arrays  are  usually 
assumed  to  represent  a  snapshot  in  time,  wiUi  only  alti¬ 
tude  differences  influencing  the  radiance  distribution.  As 
is  seen,  the  sun  obscuration  is  substantially  different  in  the 
array  representations,  and  thus  a  more  dUfficult  con^puta- 
tional  procedure  would  be  required.  The  data  in  Bg.  5-2a 
and  5-2b  however  represent  more  consistent  conditions 
with  respect  to  solar  obscuration  and  could  be  handled 
readily  in  the  traditional  computational  scheme,  DunCley 
era/.  (1976). 

Flight  C-410  produced  ten  retrievable  arrays,  six  of 
which  are  illustrated  in  Fig.  5-3.  This  series  of  plots 
represents  the  second  half  of  the  total  data  flight  and  is 
chosen  to  illustrate  the  influence  of  changing  doud  condi¬ 
tions.  At  the  be^nniitg  of  the  sequence,  Le.  during  the 
400m  AGL  runs,  one  notes  the  unobscured  sun  in  both 
spectral  bands,  plus  the  chanced  terrain  reflectance  at  the 
Filter  5  (765nm)  wavelength  (Ref.  Fig.  2*8).  At  the  inter¬ 
mediate  altitude,  the  sun  becomes  obscured  at  some  time 
between  the  Filter  4  and  the  Filter  5  setpiences,  and  by 
the  time  the  top  altitude  runs  are  made,  the  sky  is  fully 
overcast  above  300(kn  AGL.  In  all  cases  the  relatively 
high  terrain  radiances  at  765nm  are  readily  apparent. 

Flight  C-468  represents  a  dual  (2+4)  mission  in 
which  all  but  one  of  the  anticipated  sixteen  data  arrays 
were  successfully  retrieved.  One  notes  the  unobscured 
sun  occurring  in  each  display  and  the  dramatic  change  in 
the  relationship  between  sky  and  ;rrain  radiance  levels  as 
the  altitude  of  measurement  increases  allowing  the  path 
radiance  and  terrain  reflectance  properties  to  assert  their 
influence  in  the  measurements  of  the  lower  hemisphere. 


FLIGHT  €-392  -  1  NOVEMBER  1976 

DESCRIPTION  OF  FLIGHT 

AND  WEATHER  CHARACTERISTICS 
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Flltfrt  DescrtpllMi.  Flight  C-392  was  an  afternoon 
fki^t.  Multiple  doud  layers  varied  from  broken  to  over¬ 
cast.  ITie  approximate  northeast  to  southwest  Meppen 
track  was  located  between  Oldenburg  and  Lathen  in 
northweMm  Owimny.  IVpical  terrain  femtm  were 


heavily  cultivated  low  lying  flat  farmlands  interspersed 
with  occasional  dark  woods  and  small  towns. 

la*I1iglit  Nelea.  The  in-flight  observer  reported 
scattered  douds  at  1500  metm  (5000  feet),  7/8  dtocu- 
mulus  at  5400  meters  (18,000  feet),  7/8  dmis  at  6000 
meters  (20,000  feet).  By  1245  GMT  the  lower  douds 
were  scattered  to  broken  at  2100  meters  (7000  feet),  the 
altocumulus  was  unchanged  and  the  anus  had  incressed 
to  overcast.  Conditions  worsened  at  the  east  end  with  the 
lower  broken  layer  having  tops  of  2400  meters  (8000  feet) 
and  increasing  to  overcast  by  1312  GMT.  Moderate  haze 
was  observed  throughout  the  track  and  the  period. 

Local  Weatim  Nates.  Twente,  west  of  the  track, 
reported  1/8  stratocumulus  at  1500  meters  (5000  feet) 
and  7/8  altocumulus  at  2400  meters  (8000  feet)  at  KXX) 
and  1  too  GMT.  The  stratocumulus  lowered  to  750  meters 
(2500  feet)  by  1200  GMT  and  the  altocumulus  raised  to 
3600  meters  (12,000  feet).  By  1229  GMT  the  stratocu¬ 
mulus  increased  to  6/8  coverage  at  630  meters  (2100  feet) 
aiKl  gradually  lowered  to  540  meters  (1800  feet)  by  \AQO 
GMT.  The  altocumulus  layer  increased  to  overcast  and 
gradually  lowered  to  2700  meters  (9000  feet)  through  the 
afternoon.  Visibility  was  6  to  8  kilometers  in  haze. 

Lingen,  south  of  the  track,  reported  5/8  stratocu¬ 
mulus  at  1200  meters  (4000  feet)  and  6/8  altocumulus  at 
2700  meters  (9000  feet)  with  visibility  4.2  kilometers  in 
light  fog  at  1100  GMT.  Conditions  improved  to  3/8 
cumulus  and  stratocumulus  at  450  meters  (1500  feet),  S/8 
altocumulus  at  3000  meters  (10,000  feet),  visibility  7  to  8 
kilometers  in  haze.  At  1500  there  was  5/8  stratocu¬ 
mulus  at  750  meters  (2500  feet),  overcast  altocumulus  at 
2700  meters  (9000  feet)  and  visibility  10  kilometers. 

Oldenburg  and  Ahthom,  east  of  the  track,  indude 
ceiling  and  visibility  data  in  their  reports  but  no  doud 
types  or  amounts.  At  Oldenburg  the  ceilings  were 
reported  from  1200  to  1500  meters  (4000  to  5000  feet) 
and  visibility  5  to  7  kilometers  with  haze.  Ahihom 
reported  ceilings  1500  to  2700  meters  (5000  to  9000  feet) 
and  visibility  6  to  7  kilometers  in  haze. 

The  radiosonde  station  at  Rheine/Waldhugel  was 
approximately  81  kilometers  south  of  the  flight  track 
center  point,  but  no  data  were  taken  on  this  date. 

SyMgtle  RMiarka.  The  surface  chart  for  1200 
GMT  had  a  weak  ridge  with  its  axis  through  eastern  Ger¬ 
many.  From  a  960-milttbar  low  centered  south  of  fodamf 
an  oedusion  extended  east  and  south-southeast  through 
the  western  part  of  the  North  Sea  then  as  a  cold  front 
south  and  south-southwest  through  western  Ftanoe  and 
northwestern  Spam  and  Porti^  into  the  Aflamic.  At  SOO 
miRibtrs  there  was  a  low  over  western  Plfiiand.  A  weak 
gradient  prevailed  ovet  western  Europe  with  tight  to 
moderate  westerly  winds.  The  airmsss  was  stable  mari¬ 
time  polar. 
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Fig.  5-2  (con*t).  Sky  and  terrain  radiances  -  Meppen,  Germany.  Flight  C-392,  A  —  557nm  and  765nm,  0^  —  70^. 


FLIGHT  C.410  -  4  JULY  1977 
DESCRIPTION  OF  FLIGHT 
&  WEATHER  CHARACTERISTICS 
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Flicht  DescripHon.  Flight  C-410  was  a  midday 
flight,  spanning  focal  apparent  noon.  There  were  multiple 
cloud  layers  varying  from  scattered  cumulus  to  overcast 
cirrus.  The  approximate  east-west  Bniz  track  was  centered 
south  of  Rennes  in  northwestern  France.  Typical  terrain 
features  were  green  and  brown  fields  interspersed  with 
small  towns. 

In-Fligirt  Netcs.  The  in-flight  observer  reported 
broken  to  overcast  altostratus  clouds  at  the  beginning  of 


the  flight  with  scattered  cumulus  forming  at  1200  meters 
(4000  feet)  after  1153  GMT.  Tops  of  the  cumulus  were 
1500  meters  (5000  feet)  and  built  to  1800  meters  (6000 
feet)  before  the  end  of  the  flight.  Moderate  haze  was 
reported  at  all  altitudes. 

Local  Weather  Notes.  Data  from  Rennes,  7  kilom¬ 
eters  north  of  the  center  of  the  flight  track,  show  1/8  to 
3/8  cumulus  at  altitudes  varying  from  900  to  1380  meters 
(3000  to  4600  feet),  3/8  to  5/8  altostratus  at  3000  meters 
(10,000  feet)  and  7/8  cimis  at  6000  meters  (20,000  feet). 
Visibility  of  20.0  kilometers  decreased  to  U.2  kilometers 
by  the  end  of  the  flight. 

Nantes/Chateau  Bougon,  78  kilometers  south  of  the 
track  center  point,  reported  1/8  cumulus  increasing  to  4/8 
at  attitudes  varying  from  1200  to  1380  meters  (4000  to 
4600  feet)  and  3/8  to  7/8  altocumulus  at  3000  meters 
(10,000  feet).  Visibility  varied  from  25.0  to  11.2  kilome¬ 
ters. 

The  radiosonde  station  at  Brest  was  208  kilometers 
westnorthwest  of  the  flight  track  center  and  located  in  a 
prevailing  airflow  that  was  parallel  with  the  track. 

Synoptic  Remarks.  The  surface  charts  show  a  1024 
millibar  high  east  of  Denmark.  Most  of  France  was  in  a 
col  with  pressures  about  1014  millibars.  At  500  millibars 
there  was  a  low  pressure  trough  over  southern  France. 
The  flow  was  easterly  and  the  air  mass  was  modified 
unstable  maritime  polar. 
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FLIGHT  C-4M  -  21  AUGUST  ir78 
DESCRIPTION  OF  FLIGHT 
&  WEATHER  CHARACTERISTICS 
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Flight  Description.  Right  C'468  was  a  midday 
flight  spanning  local  apparent  noon  with  take  off  at  0847 
and  landing  at  1320  GMT.  There  were  scattered  high  thin 
cirrus  clouds.  The  approximate  northeast  to  southwest 
Meppen  track  was  located  between  Oldenburg  and  Lathen 
in  northwestern  Germany.  Typical  terrain  features  were 
green  and  brown  fields  interspersed  with  occasional  dark 
woods  and  small  towns. 

In-Flight  Notes.  The  in-flight  observer  noted  that 
there  was  very  muted  horizontal  banding  in  haze  and 
cirrus  above  -  gradual  clearing  of  haze  left  relatively  clear 
cirrus  layers  against  blue  sky.  A  weak  frontal  passage  last 
night  has  left  cirrus.  A  new  frontal  system  to  the 
northwest  is  causing  more  cirrus  to  move  in  from  the 
northwest  slowly.  At  0920  GMT  at  300  meters  (1000 
feet)  there  was  scattered  cirrus  at  6000  meters  (20,000 
feet)  with  8  kilometers  (5  miles)  slant  range  down  sun. 
On  the  climb  from  1350  to  3000  meters  (45<X)  to  10,000 
feel)  the  low  altitude  haze  top  was  1860  meters  (6200 
feet)  with  a  thin  clearer  layer  and  lighter  haze  above.  Dis¬ 
tinct  haze  top  of  all  haze  at  2400  meters  (8000  feet). 
Areas  of  clouds  are  visible  to  the  north  and  south.  At 
1005  GMT  at  3(XX)  meters  (10,000  feet)  scattered  cirrus  at 
6000  meters  (20,000  feet)  with  scattered  cumulus  far  to 
north,  light  haze  with  slant  range  16  kilometers  (10 
miles).  Cirrus  appears  to  come  down  to  3600  to  4500 
meters  (12,000  to  15,000  feet)  in  places.  On  the  climb  to 
6000  meters  (20,000  feet)  some  cirrus  trails  were  heavy  to 
below  4500  meters  (15,000  feet);  cirrus  are  in  bands  with 
20  to  30  nautical  mile  separation.  Other  cirrus  extend  up 
to  10500  meters  (35,000  feet);  basic  layer  is  at  6000 
meters  (20,000  feet).  At  1130  GMT  at  1350  meten  (4500 
feet)  there  was  scattered  altostratus  at  1500  meters  (5000 
feet)  and  scattered  cirrus  at  9(XX)  meters  (30,000  feet) 


with  moderate  to  heavy  haze  and  slant  ranfs  6.4  to  t 
kilometers  (4  to  5  miles).  At  1145  OMT  on  the  dimb  to 
3000  meters  (10,000  feet)  it  was  noUna  that  cumulus  were 
beginning.  Gnus  was  present  during  the  first  half  with  a 
moderate  wall  to  the  south.  A  new  band  or  drrus  was 
moving  in  from  the  north  at  a  higher  altitude  than  previ¬ 
ously.  The  cumulus  tended  to  blend  into  bright  haze  in 
the  up  sun  direction.  On  the  final  descent  the  cumulus 
field  was  developing  further  along  and  to  the  south  of  the 
west  end  of  the  track.  Cumulus  tops  were  2700  meters 
(9000  feet)  and  the  haze  top  was  now  closer  to  2100 
meters  (7000  feet).  Visibility  at  the  top  of  the  haze  layer 
was  very  poor.  Visibility  below  2700  meters  (3000  feet) 
seems  lower  than  previous. 

Local  Weather  Notes.  Bremerhaven,  87.1  kilome¬ 
ters  northeast  of  track  center,  reported  scattered  high  thin 
cirrus  clouds  becoming  broken  at  6900  meters  (23,000 
feet)  by  1300  GMT.  Visibility  of  5  kilometers  in  haze 
improved  to  9  kilometers  by  1300  GMT. 

Emden,  45.7  kilometers  northwest  of  the  track 
center,  recorded  1/8  to  2/8  high  thin  cirrus  becoming  S/8 
broken  at  6900  meters  (23,000  feet)  at  1200  GMT.  Visi¬ 
bility  of  4.5  kilometers  in  haze  improved  to  9  kilometers 
by  1000  GMT  and  decreased  to  8  kilometers  at  1200 
GMT. 

Bremen,  79.3  kilometers  east  of  the  track  center, 
reported  3/8  decreasing  to  1/8  high  thin  cirrus  with  visibil¬ 
ity  8  to  15  kilometers. 

The  radiosonde  station  at  Bergen  was  149  kilometers 
east  and  downstream  from  the  track.  The  raob  showed  a 
warming  trend  below  400  millibars  and  cooler  above  than 
the  sounding  taken  six  hours  earlier. 

Synoptic  Remarks.  The  surface  chart  for  0000 
GMT  showed  an  occluded  front  that  extended  from  north 
of  the  Arctic  Circle  southeast  and  south  to  near  Oslo  then 
as  a  cold  front  southsouthwest  through  the  North  Sea  to 
Belgium  and  off  the  west  coast  of  Morocco  with  a  wave 
near  Madrid.  The  front  was  approaching  the  flight  track 
and  there  was  southwesterly  flow  over  the  track.  At  1200 
GMT  the  cold  front  had  passed  over  the  track  and  was  on 
a  line  from  southern  Norway,  Bremen,  Offenbach,  Avig¬ 
non  to  Fes.  There  was  a  weak  cell  of  high  pressure  cen¬ 
tered  two  degrees  west  of  Brest.  Flow  over  the  track  was 
west-southwesterly.  At  500  millibars  at  0000  GMT  there 
was  weak  ridging  from  northern  Italy  to  northern  Scandi¬ 
navia  and  a  dosed  low  just  southeast  of  Iceland.  At  1200 
GMT  there  was  ridging  from  Spain  northeastward  to 
Latvia.  The  track  had  weak  northwesterly  flow  at  this 
level.  The  air  mass  was  modified  maritime  polar.  The 
satellite  maps  for  1259  and  1309  OMT  showed  thin  douds 
over  all  of  western  Europe  and  a  frontal  system  approach¬ 
ing  Great  Britain. 
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ZENTTH  ANGLE  (Defraes)  ZENIIH  ANGLE  (Depwi) 


Fig.  S-9  (cM’t).  Sky  and  terrain  Radianoes  -  Meppen,  Germany.  Flitfit  C-4M,  altitude  —  140lkn  AGL. 


ZENITH  ANGLE  (Degreea)  ZENITH  ANGLE  (Degreea) 


F1g.S-d.  Sky  and  terrain  Radianeea  •  Meppen,  Germany.  Right  C-d«,  altitude  ~  3  lOOtn  AOL. 
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Ftf.  5-7.  Sky  and  terrain  Radiances  -  Meppen.  Germany.  Flight  C-468,  altitude  —  6000m  AOL. 


6.  SUMMARY 

A  data  set  which  contains  nearly  500  arrays 
representing  measurements  of  sky  and  terrain  radiance 
values  is  described  and  illustrated.  The  measurements 
were  made  using  radiometer  systems  mounted  on  a  C-130 
aircraft  during  a  series  of  European  flights  associated  with 
the  NATO  program  OPAQUE,  Fenn  (1978)  and  Johnson 
etai  (1979). 

Radiance  measurements  throughout  the  4ir  field 
surrounding  the  aircraft  were  made  in  four  spectral  bands 
having  mean  wavelengths  of  478,  557,  664  and  765nm. 
Spectral  sets  of  measurements  were  made  at  eadi  of 
several  altitudes  between  the  surface  and  anmximately 
6  km  AGL. 

Evaluation  of  the  data  illustrates  good  orientation 
and  radiometric  accuracies,  except  for  stray  light 
influences  as  the  radiometer  field  of  view  swept  near  and 


through  the  solar  disc.  Corrective  procedures  are  outlined 
and  their  application  is  undergoing  (further  refinement. 

These  data  are  appropriate  for  use  in  the  devdop- 
ment  of  operationally  jaaiefvl  predictive  models  as  illus¬ 
trated  in  Hering  (1981).  The  radianoe  data  in  coofunction 
with  their  companion  scattering  coefBctent  data,  Johnson 
and  Gordon  (1980)  are  readily  applicable  to  the  determi¬ 
nation  of  slant  path  contrast  transmittances.  atmospheric 
optical  depths,  aerosol  directional  scattering  characteristics, 
fhx  divergences  and  their  attendant  determinations  of  tur¬ 
bid  atmos|4ieric  single  scattering  albedos. 
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APPENDIX  A:  MHcmiogical  GiMsarr  A  AMgvHatitM 


SKY  AND  CEILING 

Sky  cover  symbols  are  in  ascending  order.  Figures  prO’ 
ceding  symbols  are  heights  in  hundreds  of  feet  above 
station.  Sky  cover  symbols  are: 

O  Clear;  less  than  0.1  sky  cover 
O  Scattered:  0.1  to  less  than  0.6  sky  cover 
dP  Broken:  0.6  to  0.9  sky  cover 
^  Overcast:  more  than  0.9  sky  cover 
-  Thin  (when  prefixed);  light  (when  suffixed) 

--  Very  light  (when  suffixed) 

•X  Partial  obscuration:  0.1  to  less  than  1.0  sky  hidden 
by  precipitation  or  obstruction  to  vision  (bases 
at  surface) 

X  Obscuration:  1.0  sky  hidden  by  precipitation  or 
obstruction  to  vision  (bases  at  surface) 


Letter  preceding  height  of  layer  identifies  ceiling  layer 
and  indicates  how  ceiling  height  was  obtained.  Thus: 

A  Aircraft 

B  Balloon  (pilot  or  ceiling) 

D  Estimated  height  of  cirriform  clouds  on  basis  of 
persistency 

E  Estimated  height  of  noncirriform  clouds 
M  Measured 

R  Radiosonde  balloon  or  radar 
U  Haight  of  cirriform  ceiling  layer  unknown 

V  Immediately  following  numerical  value  indicates  a 
varying  ceiling  (also  used  witft  varying  visibility) 

Mf  Indefinite,  sky  obscured  by  surface  base  phenomenon. 
a.g.  fog.  blowing  dust,  snow 


RELATIVE  HUMIDITY  (RH| 

R^Mitad  in  peroant  and  oomputad  horn  temperature  and 
dmi^oint. 


VISISILITY  (VV) 

Reported  in  kilometers. 


WEATHER  AND  ORSTRUCTION 
TO  VISION  SYMROLS 


A 

Hail 

IF 

loa  fog 

AP 

Small  hail 

K 

Smoke 

BO 

Blowing  dust 

L 

Drizzle 

BN 

Blowing  sand 

R 

Rein 

BS 

Blowing  snow 

m 

Alin  showers 

0 

Dust 

s 

Snow 

E 

Sleet 

S6 

Snow  grains 

EW 

Slaat  showers 

SP 

Snow  pellets 

F 

Fog 

SR 

Snow  showers 

GF 

Ground  fog 

T 

Thunderstorms 

H 

Haze 

ZL 

Freezing  drizzle 

1C 

Ice  crystals 

ZR 

Freezing  rain 

CLOUD  ARRREVIATIONS 


Ac 

Altocumulus 

Cs 

Cirrostratus 

At 

Altostratus 

Cu 

Cumulus 

Cb 

Cumulonimbus 

Ns 

Nimboatratua 

Cc 

Cirrocumulus 

Sc 

Siratocumulua 

Ci 

Cirrus 

91 

Stratus 

WIND 

Direction  in  ten's  of  degrees  from  true  north,  speed  in 
meters  per  second  (mps).  A  "0000"  indicates  oelm.  A 
"G"  indicetet  gusty.  A  "Q"  indicates  squall.  Psak 
•paad  of  gusts.  uAan  reported,  follows  G  or  Q.  Ths 
contraction  NSHFT  in  remarks  followad  by  thaa  group 
(GMT)  indicatea  wind  shift  and  its  time  d  oocurranca. 

Dcaa^:  010B  is  010  degraas.  •  amo. 

aaOTGII  iaaoOdigraat.  Tmpa.paak 
apaad  in  giiata  af  It  mpt. 


•Ms  Data  miaaing  in  original  aouroa. 


APPENDIX  B  (E«r.  APGL-TB*il-flM) 


SUMMARY  OF  INPUT  DATA, 
COMPUTATION  STEPS  AND 
RESULTANT  DATA 


B.1  General  Inpnt  Deta 


B.3.1  Compute  path  tfinsmittanoe  distribution  for  each 
sub-layer. 

B.3.9  Compute  sky/terrain  radianoe  distribution  at 

sub-layer  boundaries. 

B.3.10  Compute  slant  path  contrast  transmittance  of 
selected  paths. 


B.1.1  Number  of  atmospheric  layers 

B.1.2  Representative  wavelength  (41m) 

B.  1 . 3  Extraterrestrial  solar  (lunar)  irradiance 

B.1.4  Rayleigh  optical  depth  for  top  layer  for 
\  -  .SS^m  (unless  otherwise  specified,  set  equal 
to  .0258  for  layer  base  altitude  of  10  km). 

B.l.S  Solar  zenith  angles  (deg.) 

B.1.6  Zenith  viewing  angles  (deg.) 

B.1.7  Scattering  angles  between  viewing  path  and  sun 
(deg.,  max.  10) 

B.1. 8  Average  reflectance  of  underlying  surface  and 

reflectance  of  background  in  the  vicinity  of  target 
if  different  than  the  average  reflectance. 

B.2  Input  Data  Each  Layer 

B.2.1  Base  altitude  (km) 

B.2. 2  Scattering  mixing  ratio 

B.2.3  Single  scattering  albedo 


B.4  Output  Tables  aud  Duta 


B.4.1  Single  scattering  phase  function,  as  calcu¬ 

lated  or  specified  for  each  layer. 

B.4.2  Pith  radiance,  L*(x,d,^),  for  each  sub-layer  and 
eadt  solar  zenith  angle,  P,. 

B.4.3  Path  transmittance  rr(z,P)  for  each  sub-layer. 

B.4.4  Sky  radiance,  £r(r,P,/9),  for  0*<P<85*  for 
selected  levels  (max  5)  for  each  P,. 

B.4.S  Sky  radiance,  Lr(z,P,/9)  for  95*<P<180*  for 
selected  levels  (max  5)  for  eadi  P,. 

B.4.6  Sky  plus  terrain  radianoe  Lr(z,P,^)  for 
9S”<P<180*for  selected  levels  (max  S)  for  each 
Ps 

B.4.7  Contrast  transmittance,  Tcfz ,P,/9)  between  obfect 
and  sensor  altitudes  (max.  5  altitude  intervals) 
for  0*4P^85*^and  for  each  P,. 

B.4.8  Contrast  transmittance  TAz.Bfi)  between  obfect 
and  sensor  altitudes  (max.  S)  for  95^<P<I80* 
and  for  each  P,. 


B.3  Computation  Steps 


B.3.1  Compute  Rayleigh  atmosphere  optical  thickness 
for  each  layer 

B.3. 2  Compute  total  optical  thickness  for  each  layer. 

B.3.3  Compute  optical  depth  for  the  base  of  each  layer. 

B.3.4  Compute  single  scattering  phase  function  and 
asymmetry  factors  for  each  layer.  (These  param¬ 
eters  may  be  specified  if  desired  in  lieu  of  the  use 
of  model  algorithms). 

B.3.S  Transform  optical  depth,  single  scattering  albedo, 
and  asymmetry  factor  in  accordance  with  delta- 
Eddington  approximation. 

B.3.6  Compute  path  function  distribution  for  the  base 
of  each  sub-layer  (each  layer  is  divided  equally 
into  2  sub-layers). 

B.3.7  Compute  path  radiance  distribution  for  each  sub¬ 
layer. 


B.5  Additieual  Output  Fraoi 
The  Supplemental  Eddlugten 
Computer  Program 

The  standard  Eddington  computer  program  is 
employed  as  an  integrated  part  of  the  model  calculations 
of  contrast  transmittance.  The  following  output  from  the 
program  as  modified  by  the  ddu-Eddington  transfoima- 
tion  of  input  parameters  is  available  in  tabular  form  for 
the  base  of  each  sub-layer  and  for  each  solar  zenith  angle: 

optical  depths,  r  and  r' 

components  of  Eddington  diffuse  radianoe,  L£>(z) 
and  Ld'(z) 

total  scalar  irradiance 
downward  diffuse  irradiance 
upward  diffuse  irradiance 
solar  irradiance 
total  downward  irradiance 

s 
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